Four-qubit entangled symmetric states with positive partial transpositions 
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We solve the open question of the existence of four-qubit entangled symmetric states with positive 
partial transpositions (PPT states). We reach this goal with two different approaches. First, we 
propose a half-analytical-half-numerical method that allows to construct multipartite PPT entangled 
symmetric states (PPTESS) from the qubit-qudit PPT entangled states. Second, we adapt the 
algorithm allowing to search for extremal elements in the convex set of bipartite PPT states [J. M. 
Leinaas, J. Myrheim, and E. Ovrum, Phys. Rev. A 76, 034304 (2007)] to the multipartite scenario. 
With its aid we search for extremal four-qubit PPTESS and show that generically they have ranks 
(5,7,8). Finally, we provide an exhaustive characterization of these states with respect to their 
separability properties. 



Introduction.- Entanglement has become an impor- 
tant notion in modern physics [1]. This striking feature 
of composite physical systems not only fundamentally 
distinguishes classical and quantum theories, but it has 
also developed into a key resource for various applica- 
tions. For instance, it allows for quantum teleportation 
[2], quantum cryptography [3], and is a prerequisite for 
another important resource in quantum information the- 
ory (QIT)-nonlocal correlations [4]. Deciding, then, if a 
given quantum state is entangled (i.e., if it is not a mix- 
ture of products of states representing individual subsys- 
tems [5] ) has become one of the most important problems 
(the so-called separability problem) in QIT and, even if 
simple to formulate, it is one of the hardest to solve [6]. 

Due to the recent achievements in experimental im- 
plementations of various many-body states such as, for 
instance, the four-qubit bound entangled Smolin state 
[7] , the six-qubit Dicke state states [8] , or the eight-qubit 
Greenbcrgcr- Horne-Zeilinger (GHZ) state [9], the sepa- 
rability problem in quantum systems consisting of more 
than two constituents has gained importance. Here, the 
problem becomes even more complicated because one 
wants to answer not only the simple question of whether 
a particular state is entangled, but also what type of en- 
tanglement it has (see Rcf. [10]). Various approaches 
have been proposed to detect and characterize entangle- 
ment in such systems (see, e.g., Refs. [11-14] and a recent 
review [15]). 

With this paper we fit into the above line of research 
and start a general program of characterization of entan- 
glement properties and correlations of an important class 
of multipartite states — the so-called symmetric states 1 . 
These states have already been investigated (see, e.g., 



Refs. [16, 18-21]). More attention, however, has been 
devoted to pure states, while entanglement properties of 
mixed states are mostly unstudied. In particular, it re- 
mains uncertain if there exist four-qubit entangled sym- 
metric states with all partial transpositions positive or, 
in other words, whether the separability condition based 
on partial transposition is necessary and sufficient in this 
case. It is known that PPT symmetric states of three- 
qubits are all separable [16], and existence of such states 
of five and six qubits has recently been reported [19, 20]. 
The main aim of the paper is to fill in this gap by showing, 
contrary to common belief, that there exist four-qubit 
PPTESSs. Then, we thoroughly study the entanglement 
properties of four-qubit PPT symmetric states. 

Preliminaries and general entanglement properties of 
four-qubit symmetric states.- Let us start from a couple 
of definitions that will frequently be used throughout the 
paper. Consider a product Hilbert space Hj.n = (<D d )® N 
and a convex set D of iV-partite states p acting on Hd,N- 
By r(p), K(p), and R(p) we will be denoting the rank, 
kernel, and range of p. Also, the notations (/o, . . . , fd-i) 
and |ei, . . . , ejv) will be used to denote a vector |/) G (D d 
and a pure product vector from Hd,N, respectively. 

We say that p is fully separable iff it can be written in 
the following form [5]: 
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where A\ , . . . , An denote parties and p l A . are density ma- 
trices representing respective subsystems. 

Then, splitting the parties A\ , . . . , An into two disjoint 
subsets S and S, we say that p is PPT with respect to the 
bipartition S\S if and only if p Ts > (note that the par- 
tial transpositions with respect to S and S are equivalent 



1 In other words, states describing bosonic systems consisting of 
finite number of two-dimensional subsystems. Here we assume 
the usual definition of separability, with respect to the Hilbert 
space being a product of single-particle Hilbert spaces. Other 



definitions are also considered in such systems (see, e.g.. 
[17] and references therein). 
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under the full transposition). States with this property 
make a convex subset D$ of D. An intersection of these 
subsets for all S together with D defines the set, denoted 
Z?ppt, of states that have all partial transposes positive, 
further called PPT states. An important class of PPT en- 
tangled states are the so-called edge states [22, 23]. We 
say that a given PPT state p is edge if there is no product 
vector |ei, . . . , ejv) € Hd,N such that |d, ■ • ■ , cn) G R{p) 
and its partial conjugation with respect to S is in R(p Ts ) 
for all S. Edge states are an important tool in the char- 
acterization of PPT entangled states because every PPT 
state can be decomposed into a convex combination of 
an edge state and a separable state, or, more precisely, 
edges states are those from which no separable states can 
be subtracted without loosing the PPT property. Thus, 
they are crucial for the full characterization of entangle- 
ment in PPT states. More attention has been devoted 
to edge states in biparite and three-partite systems (see 
Refs. [23, 24]), while little is known about them in TV- 
partite systems. A very convenient way to classify and 
characterize edge states is to use their ranks together 
with the rank of all the relevant partial transposes, i.e., 
(r(p), r(p TA i ),..., r(p TA N ), r (p TA i A 2 ), . . .). 

We can now pass to the TV-qubit symmetric states. 
Let us focus on i?2,iv an d denote by Sn and Pn the 
symmetric subspacc of H 2i n and a projector onto Sn- 
Recall, that Sn is spanned by the unnormalized vectors 
\Ef) = |{0, TV — *}, {1, *}>, where |{0,TV- »},{l,t}) is a 
symmetric vector consisting of i ones and TV— i zeros. For 
further benefits, let us notice that dimiSjv = TV + 1 and 
hence Sn is isomorphic to <C N+1 . We then call a density 
matrix p acting on H 2t N symmetric iff R(p) C Sn- 

In the case of symmetric states the number of rele- 
vant partial transpositions defining the set of PPT states 
can be significantly reduced. Clearly, positivity of a par- 
tial transposition with respect to a particular S implies 
positivity of partial transposition with respect to all sub- 
sets S with the same number of parties. Together with 
the equivalence of some of partial transpositions under 
the full transposition, this results in only [TV/2J of rele- 
vant partial transpositions. For concreteness, we choose 
them to be Ta i; Ta 1 a 2 > Ta 1 a 2 a 3 , etc. In the particular 
case of TV = 4 there are only two of them, which, break- 
ing the general notation, we will be denoting by Ta and 
Tab- Consequently, the set of four-qubit PPT symmetric 
states -Dp 5 ^ is an intersection of three sets D sym , D s J m , 
and D s Jg. Accordingly, one has in this case only three 
relevant ranks (r(p) , r(p TA ) , r(p TAB )) , which, for the sake 
of simplicity, we call three-rank of p and denote as r(p). 
Notice also that the fact that p is symmetric imposes non- 
trivial bounds on r(p), r(p T -*), and r{p TAB ). First of all, 
dimS^ = 5 implies r(p) < 5. Then, since S3 is isomorphic 
to (D 4 , while S 2 to (D 3 , r(p TA ) < 8 and r(p TAB ) < 9. 

Passing to the separability properties of PPT qubit 
symmetric states, it is known that for TV = 2 and TV = 3 
all are separable. While the first case directly follows 
from the results of [25], for the second one, one uses the 
fact that 52 is isomorphic to (D 3 and therefore p can be 



seen as a PPT qubit-qutrit state. Again, the results of 
Ref. [25] apply here. Finally, it was shown in Ref. [16] 
that if p is a symmetric TV-qubit state and r(p) < TV, 
then it takes the form (1). The first nontrivial, and so 
far unsolved case appears for TV = 4. All PPT four-qubit 
symmetric states with r(p) < 4 are separable, however, it 
has not been known whether the same holds for r(p) = 5. 
In other words, it remains uncertain whether there are 
no PPT entangled four-qubit symmetric states and the 
partial transposition provides a necessary and sufficient 
criterion in this case. Our main aim is to show that 
this is not the case and there do exist examples of PPT 
entangled states supported on 54. 

Before getting to the construction, let us first dis- 
cuss separability properties of four-qubit PPT symmet- 
ric states and single out all instances with respect to the 
three-rank when there are edge states. All the theorems 
proven below are left with sketches of proofs, while their 
detailed versions may be found in appendix A and Ref. 
[26]. 

Together with the already mentioned results of Ref. 
[16] we have the following theorem. 

Theorem 1. Let p be a four-qubit PPT symmetric state. 
If either r(p) < A, or r(p TA ) < 4, or r(p TAB ) < 3 then 
p is separable, while if r(p Tji ) < 6, or r{p T - AB ) < 6, then 
generic p of such ranks is separable. 

Proof. The cases of r(p) < 4 and of r(p TA ) < 4 are proven 
in Ref. [16]. The remaining ones follow from the results 
of Refs. [27, 28], which say that a PPT state p acting 
and supported on C M <S> <D N (M < TV) of rank r(p) < 
TV is separable. In the case of r(p TAB ) < 3 one treats 
P Tab as a PPT state acting on ((D 3 )® 2 . In the cases of 
r(p TA ) < 6 and of r(p TAB ) < 6 one sees p TA and p TAB 
as bipartite PPT states acting and generically supported 
on(D 2 ®(C 6 . □ 

Theorem 2. Generic four-qubit PPT symmetric states 
of a given three-rank different from (5, 7, 7) and (5, 7, 8) 
are not edge. 

Proof. Roughly speaking, for a generic symmetric PPT 
state p of a particular three-rank r(p), except for (5, 7, 7) 
and (5, 7, 8), we find a symmetric product vector |e)® 4 G 
ff 2 ,4 in the support of p such that |e*)|e)® 3 € R(p TA ), 
and |e*,e*,e,e) e R(p TAB ). Clearly, due to theorem 1 
most of the cases with respect to r (p) are already ruled 
out, (5,8,9) is trivial, and those that need to be treated 
separately are (5, 8, 7), (5, 8, 8), (5, 7, 9) (see appendix A). 

□ 

Although using the above method we cannot prove 
that generic states of ranks (5,7,7) are not edge, it is con- 
jectured to be the case. More importantly, there is an in- 
dication that they are generically separable, but this will 
be studied elsewhere [26]. Also, exploiting the methods 
below we obtained examples of edge PPTESSs of ranks 
(5,7,8); all the found examples of ranks (5,7,7) were sep- 
arable. 



Finally, let us prove that any entangled element of 
^ppt can b e decomposed in terms of at most six vec- 
tors of Schmidt rank two, i.e., entangled vectors that can 
be written as a sum of two fully product vectors. 

Theorem 3. Any entangled four-qubit symmetric state 
p can be written as 
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where K < 6, {a k ,b k ) £ (D 2 , A k ,B k £ (D. and by [ip] we 
denote a projector onto \tp). 

Proof. Applying a nonsingular transformation 
V{a,b)® 2 = (a 2 ,6 2 ) with a,b £ (D to the last two 
qubits of a four-qubit symmetric PPT entangled state 
p, one brings it to a three-qubit PPT state a acting on 
52 ® (D 2 . The latter is clearly separable and therefore 
can be written as a convex combination of at most six 
product rank-one projections [29], i.e., 
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where the form of \e k ) € 52 can determined from 
the orthogonality to K{o~) and is given by \e k ) = 
A k (a 2 k ,a k b kl b 2 k ) + B k (a 2 k ,-a k b k ,b 2 k ), while (D 2 3 \f k ) = 
|a 2 ,6 2 ) with a k ,b k ,A k ,B k £ (D. To obtain (2) and com- 
plete the proof, one applies another full rank transforma- 
tion W(a, b) (g> (a 2 , b 2 ) = (a, 6)® 3 a, b £ € to the last two 
qubits of a. □ 

Using the methods developed in Ref. [29], one can 
obtain similar decomposition in which vectors {a k ,—b k ) 
are replaced by (0,1) or (1,0) (see appendix A for the 
proof). 

Constructing four-qubit PPT entangled symmetric 
states.- We start by introducing a class of qubit-qudit 
PPT entangled states being a direct generalization of 
the 2 ® 4 PPT entangled states introduced by Horodecki 
[30] (see Ref. [31] for generalizations of 3 <E> 3 Horodecki 
states). To this end, consider the density matrices 
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|1,0>(1,0|, (4) 



where |*,) = (l/v/2)(|0,i) + |1, i + 1)) (* = 0, . . . ,d - 2). 
Then, analogously to [30], for any d > 2, we define 
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(0 < 6 < 1) (5) 



with |$ b ) = (l/V2)|0)(vT=6|0) + s/T+b\d - 1)). For 
d = 4, Eq. (5) gives the original states of Ref. [30]. 

Let us briefly characterize this class, and show that it 
maintains the separability and PPTncss properties of the 
original 2 (g> 4 Horodecki state (see appendix B for more 



details). First, one checks that for d > 2 and b £ [0, 1], 

Pib = (-"-2 ® U)pd,b{^-2 ® U^), where U is an antidiag- 
onal unitary operation consisting of unities, and there- 
fore pd,b is PPT. Second, following the argumentation 
of Ref. [30], one can show (cf. appendix B) that for 
d > 4 and b £ (0, 1), the states (5) are PPT entangled, 
and more importantly, they are edge, while for b = 
or b = 1, or d = 2, 3 separable. Third, one has that 
r(pd,b)=r(pT* b ) = d+l, 

With the aid of the states p^j, we can construct 
PPT symmetric entangled states. We will do that in 
few steps. First, we apply a full-rank transformation 
F = l d - y|0)(d - 1| with y = y/(l - b)/(l + b) to 
the second subsystem of pd,b so that the product vec- 
tors in the range of the resulting (unnormalized) state 
p' db = (1 ® F)p d , b (l <g> Ft) are given by 

(l,a)®(a d - 1 ! ...,a,l) (a € €) (6) 

(0,1)®(1,0,...,0). (7) 

The above form of the product vectors in R{p' d b ) is a key 
feature here because it allows one for a simple mapping 
of our states to many-qubit symmetric states. 

Second, to the same subsystem, we apply a nonsingular 
d' x d (d' < d) matrix 



d' d—d 
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i=0 j=0 
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(8) 



where 7* (i = 0, . . . , d — d') arc some complex parame- 
ters. By doing so, we obtain another class of PPT states 
Pd',b = (1 ® ^2)Prf. fc (l F 2 f ) which act on (D 2 (g) C d ' with 
d < d' but have additional d — d' + 1 parameters 7, . The 
transformation F 2 is chosen in such a way that it allows 
to introduce additional parameters preserving the form 
of the product vectors in the range of the resulting states. 
Precisely, the product vectors in R(p d , b ), although living 
in a smaller-dimensional Hilbcrt space, are of the form 
(6). Noticeably, since we are only using local operations, 
the resulting states are also edge. 

Third, we apply yet another local operation, denoted 
V, which maps the local vectors (a d ~ l , . . . , a, 1) to 
(1, a)®( d_1 ). Clearly, V is of full rank because in both 
vectors the same powers of a appear. By applying V to 
the second subsystem of p" d , b we simply obtain (N = d')- 
qubit symmetric states lun, which by the very construc- 
tion, have all one-particle partial transpositions positive. 

As a result the above three filters allow us to construct 
a family of many-qubit symmetric entangled states from 
the generalized Horodecki states, which are our start- 
ing point for searching for PPT entangled symmetric 
states. However, the states lon have in general nonpos- 
itive partial transpositions except for the single-partite 
ones. To overcome this, we can consider another state 
lun,\ = lun + XPn, where P/v stands for the projector 
onto the symmetric subspace Sn. Clearly there exists 
the smallest A, denoted A*, such that w;v,a* is PPT. Al- 
though it keeps the rank of the state constant, this op- 
eration, however, inevitably increases the ranks of the 



partial transpositions most probably destroying the en- 
tanglement of the resulting states. To lower them we can 
search for symmetric product vectors \e)® N G if 2, n such 
that their respective partial conjugations belong to all the 
ranges of the relevant partial transpositions R(ujJ x ) (i = 
A, AB, ...,). Every such vector can be removed from 
wjv,Ai i- e -> we consider a state ujn.x = un,\ — M-P| e )® w ! 
where P\ e \®N denotes a projector onto |e)® . By prop- 
erly choosing p and the product vector \e)® N provided 
they exist, we can lower some of the ranks of partial 
transpositions of cun,\- 

Let us now follow the above general recipe and get 
the aforementioned four-qubit symmetric PPT entangled 
states. To this end, we take p' 5 b = (1 F)p 5 ^ b (l ® Ft) 
[cf. Eq. (5)], and, following the above description, ap- 
ply the local filter F2, which is now 4x5 matrix with 
two parameters 71 and 72, and subsequently, the next 
filter V. This results in a family of four-qubit symmetric 
states £^4,6^,^ parameterized by b and 7, (i = 1,2). In 
order to get a particular example four-qubit PPTESS, 
let us put 6=1/2 and 71 = j^ 1 = l/\/2, which leads 
us to the state W4 such that uj 4 a > 0, while uj a ab ^ 0. 
To "cover" the negative eigenvalues of w 4 AB , we consider 
W4,A = w 4 + AP4. With the aid of numerics one finds 
that A* « 0.94842 is the smallest A for which uJ^ B > 0. 
However, the three-rank of W4,a, is (5, 8, 8). We can then 
lower the second rank by subtracting a product vector 
\e) m G fl(w4,A.) such that |e*)|e)® 3 G R(^ X J and 
|e*,e*,e, e) € i?(w 4 j/*). Again, exploiting numerics, we 
find that there exists such a vector |e) = (l,a), where 
a w 7+38.52091L Then, one checks that for /i « 0.64625, 
w 4.a, — M-P| e >® 4 after normalization is the expected exam- 
ple of four-qubit symmetric PPT entangled state with 
three-rank (5,7,8). Using the algorithm below one can 
check that the state is extremal, and thus also edge. It 
should be noticed that the above choice of parameters 
7j , b was made for simplicity, but other choices can also 
lead to PPT entangled states (e.g., 71 = 3/8, 72 = 11/23, 
6=1/6). 

Searching for extremal PPT entangled four-qubit sym- 
metric states.- We have just shown that four-qubit sym- 
metric PPT entangled states exist. Clearly, there must 
then exist extremal entangled elements in the corre- 
sponding set of PPT states Dp y ^. Our aim now is to 
search for such states and characterize them. For this 
purpose we adapt the algorithm for searching of extremal 
elements in the set of PPT states, originally proposed for 
bipartite systems [32] (see also Ref. [33]), to the multi- 
partite scenario. Then, we will apply it to Dpp^. 

Let us consider again the Hilbert space ifd,jv = 
(<C d )® N and the set D PPT of all PPT states acting on 
Hd,N- Let p G -Dppt and let V and V k denote projec- 
tors onto, respectively, R{g) and R{p Tk ) (fc = 1, . . . , M), 
where M denotes the number of independent partial 
transpositions (recall that some of them are equivalent 
under the full transposition). 

The state p is extremal in -Dppt iff it cannot be written 



as p = ppi + (1 —p)p2 for some pi G -Dppt and < p < 1. 
This is equivalent to say that p is not extremal iff there 
exists a density matrix a ^ p such that R(cr) C R(p) 
and R(a Tk ) C R(p Tk ) for all k. One can even relax the 
assumption of a being positive to be Hcrmitian and such 
Hermitian matrices are solutions to system of equations 



V k {h) = h (fc = 0,...,M), 



(9) 



where V = V{-)V and V k {-) = [V k {-) Tk V k ] Tk (k = 
1, . . . , M). This system is equivalent to the single equa- 
tion [Vm o . . . o V\ o V\(\i) = h. Clearly, p is a partic- 
ular solution of the system (9) and, due to the above 
statements, is extremal iff it is its only solution, which 
gives necessary and sufficient condition for extremality 
[32, 33]. This also leads to a simple necessary criterion 
for extremality. Precisely, each equation in (9) imposes 
(pN _ j r ^r fc ^j2 i mear constraints on the matrix h. The 
maximal number of conditions imposed by the system (9) 
is then ^2 k (d 2N — [r{p Tk )} 2 )- Then, a Hermitian matrix 
h has d 2N real parameters, and hence if 



M 



Y,lr(p n )?>Md 2N + l, 



(10) 



k=0 



the state p is not extremal. 

All this induces a method of searching for the extremal 
states in -Dppt [32]. Taking p G -Dppt, one solves the 
corresponding system (9). If the latter has a solution 
h 7^ p, the state p is not extremal. One then considers a 
family of matrices p(x) = (1 — xTrh)p + xh with x be- 
ing in general a real parameter. Clearly, there is x = x* 
such that P2 = p(x*) is a still PPT state, however, either 
r (/°2) = r(p) — 1 or r(p 2 k ) = r(p^) — 1 for some k. We 
can again apply the above procedure to p2, and in case 
it is not extremal get another PPT state p% with at least 
one of the ranks diminished by one. We keep applying 
this procedure until we obtain an extremal state, which 
appears in a finite number of repetitions. If the resulting 
state is pure, it is separable, otherwise it is entangled. 
Notice that in order to get a particular extremal entan- 
gled state one has to properly choose the initial state 
which basically means that it should be of higher ranks, 
as for instance the maximally mixed state, and the direc- 
tions which follow from solving Eqs. (9). 

Let us apply the above algorithm to the symmetric 
states. In this case the left-hand side of (10) has to be 
modified as p and p Tx (X = A, AB, . . .) are supported on 
Hilbert spaces of different dimensions. In particular, for 
N = 4 such analysis was done in Ref. [33] and it follows 
that states of ranks (5, 7, 9), (5, 8, 8) cannot be extremal. 
Moreover, theorems (1) and 2 imply that generic states 
of ranks (5,8,7) and (5,r(p TA ),r(p TAB )) with r(p TA ) or 
r(p TAB ) < 6 also cannot be extremal. The natural can- 
didates for extremal states have then ranks (5, 7, 7) and 
(5,7,8). 

We applied the above algorithm to four-qubit PPT 
states and all the extremal examples we found have ranks 



(5,7,8). In 30 000 runs we generated 5760 unitarily 
nonequivalent extremal entangled states and all of them 
have ranks (5, 7, 8). As an initial state we took the pro- 
jector P\ onto 5*4 (recall that the initial state has to be 
of rank five and due to theorem 1 must also have appro- 
priately high ranks of p Tx (X = A, AB)). At each step 
of the algorithm we used solutions of (9) chosen so that 
we could reach one of the three-ranks not excluded by 
theorem 2. We also got 24 240 states of ranks (5, 7, 7) in 
this way but they all are separable. 

Conclusion.- The main aim of this note was to solve 
the open question of the existence of four-qubit PPT en- 
tangled symmetric states. We have reached this goal 
by proposing a half-analytical-half-numcrical method of 
constructing of such states. The analytical part of the 
method allows one to map a class of qubit-qudit PPT 
entangled states onto many-qubit entangled symmetric 
states. Then, using already-well-established methods of 
the theory of entanglement, and with the help of numer- 
ics, we have found the desired PPT entangled states. 

We have also characterized the four-qubit PPT sym- 
metric states with respect to separability, edgeness and 
cxtrcmality properties. First, we have proven that 



generic states of a given three-rank different from (5, 7, 7) 
or (5, 7, 8) are not edge. Then, by adapting to the multi- 
partite case an algorithm allowing to search for extremal 
PPT states [32], we have sought extremal four-qubit PPT 
entangled symmetric states. All the entangled states 
found in this way have ranks (5,7,8), while those with 
ranks (5, 7, 7) encountered in this way are separable (see 
Ref. [26] for more details). 

Interestingly, all the methods presented in this paper 
can be applied to iV-qubit symmetric states. For in- 
stance, we have shown that with our method one can 
obtain five-qubit and six-qubit PPT entangled symmetric 
states, confirming the findings of Ref s. [19, 20]. Gener- 
alization of these findings to -/V-qubit symmetric Hilbert 
spaces is currently being studied and will be a subject of 
a forthcoming publication [26]. 
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Appendix A: Characterization of four-qubit 
symmetric states 

Here we recall theorems 1,2, and 3 and prove them in 
detail. 

Theorem 1. Let p be a four-qubit PPT symmetric state. 
If either r(p) < 4, or r(p TA ) < A, or r(p TAB ) < 3 then 
p is separable, while if r(p TA ) < 6, or r(p TAB ) < 6, then 
generic p is separable. 

Proof. Although the first two cases of r{p) < 4 and 
r(p TA ) < 4 were already proven in Ref. [16] let us, for 
completeness, recall their proofs. First, due to the fact 
that S3 = (D , one can always treat p as a qubit-ququart 
PPT state. Then, the results of Ref. [27] say that any 
qubit-ququart PPT state of rank r(p) < 4 is separable. 
Replacing then p by p T - A and following the same argu- 
ments, one proves the case of r(p TA ) < 4. 

In order to prove the case of r(p TAB ) < 3, one con- 
siders a state a = p TAB and exploits the fact that S2 is 
isomorphic to (D 3 . Then, a is a two-qutrit PPT state 
such that r{a) < 3, and it was shown in Ref. [28] that 
any two-qutrit PPT state of rank less or equal to three 
is separable. 

In the case of r(p TA ) < 6 let us define a = p TA and 
consider it as a bipartite state with respect to the par- 
tition B\ACD. Clearly, in such case a acts on (D 2 ® (D 6 
and is of rank at most six. Provided that it is supported 
on (D 2 ® (D 6 , which generically is the case, the results of 
Ref. [27] tell us that a is separable across B\ACD, i.e., 



-Fa 



= }jPi\ei)(ei\ B <8> \ipi}(ipi 



ACD, 



(Al) 



where, for the time being, 1^) are entangled states from 
(D 2 <E> S*2. On the other hand, the BCD subsystem of 
a is still supported on the three-qubit symmetric sub- 
space. As a result, any vector \ei)\ipi) in the decompo- 
sition (Al) must obey P 3 \ei) b\iPi)acd = \&i) b\^i) Ac d , 
where P3 is applied to BCD subsystem. This, after some 
algebra, implies that \ij)i)ACi> = \fi)A\ei)c\ei)D for some 
\fi) G <D 2 , and hence a, and accordingly p = a TA are 
fully separable. 

To prove the last case of r{p TAB ) < 6 one follows the 
same lines as before substituting p TAB for p TA . □ 

Theorem 2. Generic four-qubit PPT symmetric states 
of a given three-rank different from (5, 7, 7) and (5, 7, 8) 
are not edge. 

Proof. We will show that in all relevant cases with respect 
to the three-rank, except for (5,7,7) and (5,7,8), there 
exists a symmetric product vector |e)® 4 G R{p) such that 
|e*)|e)® 3 G R(p Ta ), and |e*,e*,e,e) G R{ P Tab ). Clearly, 
many of the three-ranks can be ruled out with the aid of 
theorem 1, and the remaining ones are (5, 7, 7), (5, 7, 8), 
(5,8,7), (5,8,8), (5,7,9), and (5,8,9). The last one is 
trivial because all symmetric product vectors |e)® 4 be- 
long to R(p) and their respective partial conjugations to 



R(p TA ) and R(p TAB ). In what follows we give a proof for 
the cases (5, 8, 7), (5, 8, 8), (5, 7, 9). 

In the case of r(p) = (5,7,9), p and p TAB are of 
full rank and therefore one has to find a product vector 
|e* , e, e, e) which is orthogonal to the only vector \fy) from 
K(p TA ). To this end, let us write |*) = |0}|V>o) + l^l^i) 
with \^i) e S3, and take |e) = (l,a) {a g <D). The or- 
thogonality condition then reads V^a) + a*W3(a) = 
with V3 and W3 denoting polynomials of degree at most 
three over the complex field. In Ref. [29], this equa- 
tion was shown to have generically (both polynomials V3 
and W3 are of degree three) at least one solution. Conse- 
quently, generic four-qubit PPT states of of ranks (5, 7, 9) 
are not edge. 

In the case of r{p) = (5,8,8), p and p TA are of full 
rank, and so one has to find a product vector \e*e*ee) 
with |e) G (D 2 orthogonal to the only vector |\&) from 
the kernel of p\ B . For this purpose, let us notice that 
p TAB = G(p*) TAB G' , where G denotes an operator swap- 
ping subsystems AB and CD. This means that |\1/) en- 
joys the same symmetry, i.e., G\^) = |^*). As a result, 
one can express it as 



l*> 



fc=i 



Xk\e* k )AB\ek)cD, 



(A2) 



where A^ G 1R and |e&) arc orthogonal symmetric two- 
qubit vectors. Exploiting the fact that \ i $l) £ K(p TAB ), 
one sees that 



K 



K 



^2 A fc (x*, y\p TAB \e* k ,e k ) = J^ X k (e k ,y\p\x, e k ) 



fe=i 



fe=i 

K 



^2\ k (e k ,x\p\e k ,y) = 

(A3) 



fc=i 



holds for any pair of qubit vectors \x) and \y) with the 
second equality stemming from the fact that p is sym- 
metric and hence p = pG. This immediately implies that 



^2Xk{e k \p\e k ) = 0, 



(A4) 



where the right-hand side is a two-qubit matrix acting 
on the CD subspace, obtained by "sandwiching" p with 
|efc)s on the first two qubits. 

On the other hand, taking into account Eq. (A2), there 
exists |e) € (D 2 such that \e*e*ee) € R(p T - AB ) iff 



J2 Xk \ e ^ ek \ 



(A5) 



In order to show that such |e) exists, assume, in contrary, 
that Eq. (A5) does not hold for any |e) G (D. Then, its 
left-hand side must have the same sign for all |e), say 



positive (as otherwise, from continuity, there would exist 
|e) for which (A5) holds). Consequently, 



e,e 



^ >>k\ek)(ek\ 



|e,e)>0 



(A6) 



for any |e) G (D 2 , which, owing to the fact that \e k ) are 
symmetric, implies that W = J2k ^k\e k ){e k \ is a two- 
qubit entanglement witness. Since all two-qubit wit- 
nesses are decomposable, we have W = P + Q Ta with 
P, Q > 0. This, when substituted to Eq. (A4), implies 
that the two conditions 



Tr[(P ® |j/)(x|)p] = 0, Tr[(Q ig) |2/X^I)P 



,Ta1 _ 



(A7) 



must be obeyed for any \x), \y) G (D 2 , contradicting the 
fact that p and p TA arc of full rank. Notice that this 
proof is general (not generic) meaning that there are no 
four-qubit symmetric edge states of ranks (5, 8, 8). 

Let us now pass to the most involving case of r{p) = 
(5, 8, 7). Here p and p TA are of full rank, while K(p TAB ) 
has dimension two. Hence, to find a product vector 
\e) m G R(p) such that |e*, e*, e, e) G R(p TAB ), one has to 
solve two equations (e*,e*,e,e\^i) = (i = 1,2), where 
l^i) are two orthogonal vectors from K(p TAB ). Exploit- 
ing again the identity p TAB = G(p*) T - AB G\ it is fairly 
easy to see that they can be written as 



l*i) 



2 

E 

/t=i 



A fc |e fc >|/ fc * 



l* 2 



2 

E 

fc=i 



A fc |/ fe )|e£). (A8) 



To see it explicitly, let us first notice that we can assume 
that one of |\frj) is of Schmidt rank two. If both of them 
are of rank three, then there exists a vector of Schmidt 
rank two in span{|^i), l 1 !^)}- On the other hand, if one 
of l^i) (i = 1,2) is of rank one, i.e., is product with 
respect to the partition AB\CD, r(p) = 4 contradicting 
the assumption that p is entangled. Assume then |\&i) is 
of rank two. Then either G|\I/|) is linearly independent 
of |*i) leading to Eq. (A8), or G|*J) = £|*i) for some 
£ G (D. In the latter case, short algebra implies that |^,) 
(i = 1,2) are not linearly independent contradicting the 
fact that they span two-dimensional kernel of p TAB . 

As a result, finding a vector |e*,e*,e,e) G R(p TAB ) is 
equivalent to solving an equation 



V(a*)W(a) + V{a*)W{a) = 0, 



(A9) 



where V, V and W, W are polynomials generically of de- 
gree two. A solution to this equation exists if and only if 
there exists z G <D such that 



and 



V(a*) =zV(a*) 



W{a) = -zW{a). 



(A10) 



(All) 



We have then brought a single equation of the fourth 
degree to two equations of the second degree. With the 



aid of the first one, we can determine a* as a function of 
z. There are clearly at most two such solutions. Putting 
them to the second equation and getting rid of the square 
root, we arrive at 



(z*) 2 Wi(z) + z*Wi(z) + Wi'(z) = 0, 



(A12) 



where W4, W±, and W'l stand for polynomials which are 
generically of fourth degree. In what follows, we will 
show that Eq. (A12) has at least one solution z = rs 
with \s\ = 1, i.e., z* = r/s. To this end, let us consider 
two cases, when s = rx and s = x/r. Putting all this to 
Eq. (A12), one gets the following equations 



( - ) W 4 (r 2 x) + -Wi(r 2 x) + W^(r 2 x) 
\x J x 



and 



— W 4 (x) 

x ' 



— Wi(x) + Wi'(x) = 0. 

x 



(A13) 



(A14) 



In the limit of r — > 00 the first equation has two roots 
s°° — > 00 (i = 1,2), while the second one four roots 
s°° -> (i = 1, 2, 3, 4). Then, in the limit of r -> 0, Eq. 
(A13) again has two roots s° — > 0, while Eq. (A14) has 
four roots s° — > 00 (i = 1,2,3,4). Notice that, after the 
above substitution, (A12) is of sixth degree in s meaning 
that it can have at most six solutions with respect to s. 
Consequently, by varying continuously r from zero to 00 
we see that at least one of the roots s,, being a continuous 
function of r, must go from zero to 00, and so there is 
a value of s such that \s\ = 1. As a result, there is at 
least one z G <D for which Eq. (A12) is fulfilled, and 
simultaneously at least one a G (D obeying (A9). □ 

Theorem 3. Let p be an entangled symmetric PPT four- 
qubit state. Then it can be written as 

p = J2lM<>k,b k ) m + B k (a kl -b k n (A15) 

fc=i 

where (a k ,b k ) G (D 2 and A k , B k are some complex coeffi- 
cients, and by [ip] we denote a projector onto \ip) . 

Proof. First, let us introduce two linear transformations 
V : (C 2 )® 2 h-> (D 2 and W : ((D 2 )® 2 i-> (C 2 )® 3 defined as 



and 



V[(a, 6)0 (a, b)} = (a 2 , b 2 ) 



W[(a,b)<E)(a 2 ,b 2 )]^(a,b) i 



(A16) 



(A17) 



respectively, with a, b being any complex numbers. Then, 
by V and W we denote maps that are defined through 
the adjoint actions of V and W, i.e., X(-) = X(-)X' 
(X = V,W). 

The key feature of the two matrices V and W is 
that WVBc\i>) = \i } ) f° r an Y three-qubit symmetric l^), 



where BC denote qubits subject to V. The same holds 
when V is applied to any pair of qubits in \ip) and followed 
by a proper application of W. Accordingly, any iV-qubit 
symmetric state p is left invariant under a proper appli- 
cation of V and W to any three-qubits. In particular, for 
a four-qubit symmetric state p, Wbc ° Vcd(p) = P- 

Let us now consider a four-qubit PPT symmetric state 
p. By applying V to the CD subsystem of p. we get 
a three-qubit state uabc = Vcd(p) acting on S2 <g> C 2 , 
where C" denotes the qubit resulting from the application 
of V. Clearly, the map V preserves positivity of partial 
transposition with respect to the first two parties, i.e., 
a T A B > o Since S2 is isomorphic to (D 3 , results of Rcf. 
[25] imply that a is separable across the cut AB\C and 
so a takes the form 



K 



a = 5> fe }(e fc |®|/ fc )<A|, 



(A18) 



fe=i 



( a k'^k) w ith ak,bk £ (D and solving 



with |efe) € (D 3 and |/fe) £ (D 2 being in general unnormal- 
ized vectors from (D 3 and (D 2 , respectively, and K < 6 
[29]. 

By the very assumption p is entangled and therefore 
r(p) = 5, which together with the fact that r(V) = 2, 
mean that the rank of a is also five. Therefore, K{a) 
consists of a single vector \<p) 6 (D 3 <E> (D 2 , which, due 
to the fact that the range of a is spanned by the vec- 
tors (a,&)® 2 <g> (a 2 ,b 2 ), takes the form \<j>) = |01) - |20). 
As a result, any product vector in Eq. (A18) has to be 
orthogonal to | 

Putting I/O 

the equation (<f>\ek, fk) z - with respect to \eu) one 
finds that it can be written as \e k ) = A k (a 2 , a k bk, b 2 ) + 
-Bfe(a 2 , — afc6fc,6 2 ) with some A k ,B k 6 (D. Putting the 
above forms of |ejt) and \fk) to Eq. (A18), one sees that 
a can be written as 

K 

a = YsKM^bkf 2 + B k (a k ,-b k )® 2 ) ® (a 2 , 6 2 )], 
fc=i 

(A19) 
where A" < 6 and [ip] denotes a projector onto \i/j). One 
completes the proof by applying W to the last two qubits 
of a. D 



Proof. The proof exploits the method developed in Ref. 
[29]. First, one notices that any p can be written as a 
sum of rank-one matrices 



K 



p = £|*0<*« 



(A21) 



where, in particular, |\&j) can be (unnormalized) eigen- 
vectors of p, and K < 6 (see the proof of theorem 3). On 
the other hand, p can always be expressed with the aid 
of the symmetric unnormalized basis {\E^ )}^ = i spanning 
54 as 



J2 P^\E*){Et\ 



(A22) 



Both decompositions (A21) and (A22) are re- 
lated via the so-called Gram system of p, i.e., 
a collection of A'-dimensional vectors |u„) = 
(1/(^|^))((* 1 |^),...,(^|^)) Ox = 1..--.5), 
giving p M „ = (u M |tv). Putting the latter to (A22) with 
explicit forms of the vectors \Vp), one recovers (A21). 

Now, by projecting the last party onto |0) we get a 
three-qubit symmetric PPT state p, which, as already 
stated, is separable. Then, according to Ref. [29], there 
exists a diagonal matrix M = diag[a^, . . . , a* K ] such that 
|u M ) = A/ M_1 |ui) (p = 1, . . . ,4). For convenience we can 
also put 1 1)5) = M 4 \vi) + \v) with \v) being some K- 
dimensional vector. Then, putting \v\) = (A\,...,A* K ) 
and \v) = (B*, . . . , B* K ), one sees that 



l*i 



E 



(K\*i 



=1 <W 



I-K 



= A i ^ar 1 |^ fe 4 ) 



BflS 



t|-"5/ 



^(1,0*)®* + ^!^), 



(A23) 



where the second equation follows from the explicit form 
of the vectors \Vff). Substituting vectors \^i) to Eq. 
(A21), one gets (A20), which completes the proof. □ 



Utilizing the normal matrix approach to the separabil- 
ity problem [29], one can prove a bit different decompo- 
sition. 

Theorem 4. Let p be an entangled symmetric PPT four- 
qubit state. Then it can be written as 

K 

p = ^[MakM)® 4 + B k (0, 1)®% (A20) 

fc=i 

where K < 6, (a,k,bk) € (D , and Ak,B k are some com- 
plex coefficients, and by [ip] we denote a projector onto 



Appendix B: Properties of the states pd,b 

Here we characterize the states (5) in more details. In 
particular, we prove that for d > 4 and b € [0, 1] they are 
PPT entangled and edge. 

Theorem 5. The states pd b are PPT for d > 2 and 
be [0,1]. 

Proof. First, notice that we can rewrite (5) in the matrix 
form as 



1 



Pd,b = 



C bBr 



(2d -1)6+1 \bB L bl d J ' 



(Bl) 



where Bjj and Bl arc d x d dimensional matrices with 
entries 1 on the upper and lower diagonal, respectively. 
Further, C is a d x d matrix given by 



/ 1 + 6 
6 



C -2 



(.) 







\/l —r 62\ 






1 + 6 ) 



(B2) 



Clearly, the partial transposition of pb d reads 



„Ta 



C bB L 
(2d -1)6+1 \bBu bl d 



(B3) 



Now, let us consider the unitary matrix U = 
antidiag[l, . . . , 1] (anti-diagonal matrix consisting of uni- 
ties). Straightforward calculations show that UBjjU^ = 
B L , UB L W = Bu, and UCW = C. Consequently, 



Pl% = (12 <8> C/)pb,d(l 2 ® C/ f ), 

meaning that p b ^ > iff pbA > 0. 



(B4) 
D 



Theorem 6. T6e states pdA are entangled for d > 4 and 
6 G (0, 1), tu/w7e separable for d = 2, 3, or 6 = 0, or 6 = 1. 

Proof. First, let us prove that for d > 4 and 6 £ (0, 1), the 
states (5) are entangled. For this purpose, it suffices to 
use the necessary criterion for separability formulated in 
Ref. [30] - the range criterion. It says that if a given den- 
sity matrix p is separable then one is able to find prod- 
uct vectors |e, /) spanning R(p) such that |e*,/) span 
R{p TA ). In what follows we show that none of the prod- 
uct vectors |e, /) in R(p) is such that |e*,/) € R(p TA ). 
All product vectors in the range of pb,d are given by 



„<J-i 



(l,a) ®{a~ - +y,a 
(0,1)0(1,0. ..,0), 



d-2 



,a,l) 



(a e <D) (B5) 
(B6) 



where y = y/(l — 6)/(l + 6). If we allow for infinite a, 
the vector (B6) may be obtained from the class (B5). 
It is also worth mentioning that the above vectors span 

R(p d ,b)- 



On the other hand, all the vectors in the range of p b A d 
are given by 



(oi, . . . ,a d -i,yai + a d ;a 2 ,a 3 , . . . ,a d ,b) 



(B7) 



with ai,...,a<j, b € (D. Consequently, a product vector 
from the first class (B5), when partially conjugated with 
respect to the first subsystem, belongs to i?(p^), i.e., 
takes the form (B7), if and only if the conditions are 
satisfied: (i) a(l - |af) = 0, (ii) a d ~ 2 = a*^" 1 + y), 
and (iii) y(y + a d ^ 1 ) = 1 + \a\ 2 . The first condition 
is satisfied if either a = 0, which contradicts the third 
condition because y ^ 1, or \(3\ 2 = 1, which contradicts 
(ii) because y ^ 0. Along the same lines one checks that 
the vector (B6) is not of the form (B7). 

In conclusion, the states pb t d are entangled for d > 4 
and6e (0,1). 

Let us finally consider the missing cases of d = 2, 3 or 
6 = 1 or 6 = 1. For d = 2 or d = 3, theorem 5 says 
that pd.fc are PPT for any 6. It is known [25] that all 
qubit-qubit and qubit-qutrit PPT states are separable. 

For 6 = it follows from Eq. (5) that p d ,o = |$o)($o|, 
which is separable, while pd.i can be written in the fol- 
lowing separable form (cf. Rcf. [30]): 

1 f 27r 
Pd,i = -7-\ dcpP(cp)®Q(<p), (B8) 

where P(tp) and Q((p) are projectors onto (l/v / 2)(l,e' ¥ ') 
and (l/^XLe-^e- 2 ^,...^-^- 1 ^), respectively. 

□ 

Theorem 7. r(p d ,b) = r{p T d A b ) =d + l. 

Proof. Direct check shows that the vectors 

|* i ) = |0,*)-|l,* + l) (i = l,...,d-2) 



|*<) = -VI + 6|00) + VI - 6|0, d- 1) + VT+ 6|11). 

(B9) 

belong to the kernel of pd,b and the subspace they span 
has dimension d — 1. On the other hand, the product 
vector from R(pd^), given in Eqs. (B5) and (B6) span 
[d + l)-dimcnsional subspace. Consequently, r(pd^b) = 
d + 1 and, since p T d \ = (1 2 ® U)p d , b (l 2 ® C/ f ), r(pJ A b ) = 
d+1. D 



